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1. Introduction. 

It is generally believed that the most profound features of the QCD- confinement and 
nontrivial 9 dependence , the chiral symmetry breaking and the absence of the ninth 
Goldstone boson are tightly connected. Let me briefly sketch these old arguments. 

We use the standard notations and the 9 vacuum is introduced to the theory in the 
following way: 

AL = 9Q (1) 
where Q is the topological charge which can be written for the 4d gluodynamics as follows: 

Q = 3^J d * xG ^ G l» = 3^2 / d ^ K ^ a = 1, 2, 3. n,v= 1, 2, 3, 4. (2) 

Here G a is the field strength tensor. It is known , the 9 term preserves renormalizability 
of the theory but is P and T odd. As it can be seen from (0) the 9 term is a full divergence 
and therefore it is equivalent to certain boundary conditions . However as is known, the 
9 is the physical parameter of the theory and the 9- dependence of physics is linked to 
the U(l) problem. • Indeed, if we believe that the resolution of the U(l) problem 

appears within the framework of these papers ,we must assume that the correlator 

K = i [d*x< OIT^-G" G° (x), -^GL G* (0)10 > (3) 
J 1 327r 2 ^ ^ 32tc 

is nonzero in pure Yang Mills theory (YM). It means that the vacuum expectation value 
(vev) of the topological density 

< \^G; v Gl u \ >= \K9, (4) 

is also nonzero . Therefore, the 9 is observable parameter. The nonzero vev < GG >^ 
does imply the CP -violation in physical transition and leads , for example ,to the mixing 
of the heavy quarkonium levels with J p = + and J p = 0~ . Only dispersion relations 
are used to translate < GG >^ into a proof of CP- violation in physical effects. ||. 

So, if we believe that U(l) problem is solved in the framework of Witten-Veneziano 
approach we automatically get K ^ (||) and therefore, the nontrivial 9 dependence. 
With introduction of the light quarks (u,d,s) the value of correlator K(^j changes Q,but 9 
parameter is still an experimentally observed quantity. 

The next link comes from the analysis of the Ward Identities (WI) which imply |4j 
that the U(l) problem is linked to the chiral symmetry breaking phenomenon. Indeed, 
the relevant WI ( in QCD with Nf flavors) takes the form: 

2 

in — Tfi f — — 

Kqcd = — <ipip> / dxi)^(x),ip^(0) (5) 

* f 

and if we have the dynamical solution of the U(l) problem (by other words,we are in po- 
sition to calculate the singlet correlator in @ to show its smallness in the chiral limit) we 

2 It is obviously that K ~ m J irrespectively to the number of light flavors (|] . 
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should explain the relation between < ipip > condensate ,responsible for the chiral sym- 
metry breaking and topological susceptibility , responsible for the solution U(l) problem. 

Besides that , in a similar way one can check that nontrivial 9 dependence in pure YM 
theory is through 9/N at large N ( in particular, < GG >~ sin(^) )H. Such a function 
can be periodic in 9 with period 2ir only if there are many vacuum states for given values 
of 9. These vacua should not be degenerated due to the vacua transitions , however the 
trace of the enlargement number of the vacuum states have to be seen in the course of 
a dynamical solution of the U(l) problem. Therefore , we have a link between the U(l) 
problem and an additional (discrete) classification of the vacuum states in the theory 
apart from the standard 9 classification. 

Another link , between U(l) problem and confinement phenomenon , can be under- 
stood from the analysis of the effective lagrangian describing the low-energy spectrum 
and dynamics of the pseudoscalar nonet in the large N limit ||,||. The most impor- 
tant assumption which has been made in the deriviation of the corresponding effective 
lagrangian was confinement. As it is known the obtaining lagrangian perfectly describes 
all properties mentioned above . In particular this lagrangian reproduces the correct 9/N 
dependence in YM theory. 

Indeed, gluodynamics can be understood as a QCD with very large quark's mass. In 
this limit effective lagrangian was founded in |J and it turns out that the number of 
vacua is of order N at iV — > oo. This fact actually is coded in the effective lagrangian 
containing the multi-branched logarithm logdet({7). In the Veneziano approach 0] the 
same fact can be seen from the formula for multiple derivation of the topological density 
Q with respect to 9 at 9 = 0. 



Therefore,the main idea of this sketch is as follows. All problems under consideration 
are tightly connected. Thus , any selfconsistent dynamical solution of one of them should 
be necessarily accompanied by the resolution of the rest problems within same approach. 

The purpose of this letter is to demonstrate that such links indeed take place within 
framework of the dynamical toron ^ approach which was discussed early in context of 
different field theories, see ref.ffl] and references therein. I would like to recall that in 
all known cases, the toron calculations give, at least, the selfconsistent results. Thus, 
one may expect that the analogous situation should takes place in the theory under 
consideration, i.e. in 4 dimensional gluodynamics with SU (2) gauge group. 

We discuss the statistical ensemble of quasiparticles which, presumably ||, describes 
the grand partition function of the 4d YM theory and which possesses strong quasipar- 
ticle interaction. It was shown that this ensemble describes the system with nontrivial 

3 We keep the term "toron", introduced in ref.Q. By this means we emphasize the fact that the 
considering solution minimizes the action and carries the topological charge Q = 1/2, i.e. it possesses all 
the characteristics ascribed to the standard toron [H . However I should note from the very beginning of 
this paper, that our solution has nothing to do with the standard toron and it is formulated in principle 
in another way than in ref. ||. The keeping of this term has a historical origin. 




Q2n-1 
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(6) 



n = 1,2... 



2 



9- dependence irrespectively to the strength of quasiparticle interaction . It should be 
noted,that the coexistence of the strong quasiparticle interaction and the nontrivial 9 
dependence is the nice feature of YM theory . 3d QED also possesses the long range 
quasiparticle interaction and confinement phenomenon flIU| . In this simple model physics 
becomes 9 independent just because of the strong quasiparticle interaction [fllfl - Crucial 
point is a very nontrivial algebraic structure of the quasiparticle interaction in YM theory, 
in comparision with analogous calculation of ref.[|lTJ in Polyakov's model. 



We will discuss the main assumptions which have been made in the description || of 
this statistical ensemble a bit later. Now I would like to cite the result of calculation of 
the topological density as a simplest application of the formulae obtained in ref. [|J: 

< GG >= i2A 4 sin(-), -vr < 9 < tt, (7) 
2 

where the A is some renormalization invariant calculable combination f\ 

The nonzero value for topological density and its 9 dependence are in agreement with 
main assumption of Witten-Veneziano approach [|TJ,@. Indeed, the formula (^) implies 
that the topological susceptibility which is nothing but derivative of the topological density 
with respect to 9, is nonzero. Besides that, this formula is in agreement with Veneziano 
expression (P) for the multiple deriviation with respect to 9. 

As it was argued above, the dynamical calculation of the topological density should 
be accompanied by the resolution of the rest problems. One such link was established 
already in the previous papers [BJ, Namelyit was demonstrated that the calculation 
of < GG > is accompanied by appearing of the additional quantum number (apart from 
the standard parameter 9) classifying a vacuum states. 

In this letter we will consider the less trivial applications of the formulae describing 
the grand partition function of YM theory , ref. |§ . Namely, we will calculate the vev of 
the Wilson loop to demonstrate the confinement in this theory. This is the main result of 
the work. 

Besides that, we find the topological susceptibility by explicit calculation ( and not 
by differentiating < GG > with respect to 9) in order to check the selfconsistency of 
the approach. In such calculation we are able to keep the nonleading term (as k — > 0) 
in the correlator for the susceptibility. It turns out ,that this term is of order (k 4 ) and 
not of order (k 2 ) as it could be expected naively. Such behavior for the correlator is the 
direct result of the strong quasiparicle interaction and might have a phenomenological 
consequences. For example, the proton "spin" problem can be reduced to the problem of 
calculation of the first moment of the topological susceptibility JT3| . 



Before we proceed to the detail consideration of the Wilson loop and topological 
susceptibility let me briefly formulate the basic assumptions of the toron approach. 

i) I allow the configurations with fractional topological charge (one half for SU(2) 
group) in the definition of the functional integral. It means that a multivalued functions 
appear in the functional integral. However, the main physical requirement is - all gauge 



4 Formula (Q) is written in Euclidean space because its calculation based on the self dual solution 
defined in this space. In contrary, the formulae (0,0) is in Minkowski space. 
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invariant values must be singlevalued . Thus , the different cuts accompany the multi- 
valued functions should be unobservable, i.e. the gauge invariant values coincide on the 
upper and on the lover edges of the cut. 

The direct consequence of the such definition of the functional integral is the appearing 
of the new quantum number classifying a vacuum states. Indeed, as soon as we allowed 
one half topological charge, the number of the classical vacuum states is increased by 
the same factor two in comparision with a standard classification, counting only integer 
winding numbers |n >. 

Of course, vacuum transitions eliminate this degeneracy. However the trace of enlarge- 
ment number of the classical vacuum states does not disappear. Vacuum states now 
classified by two numbers : < 9 < 2tt and k = 0,1. Let me repeat that origin for this 
is our main assumption that fractional charge is admitted and therefore the number of 
classical vacuum states is multiplied by a factor two. 

I have to note that the same situation takes place in the supersymmetric YM theory, 
but in this case the vacuum states are still degenerate after vacuum transitions. The 
number k in this case just numerates different vacua at the same 9. However, the gauge 
classification for the winding vacua is the same for both models (supersymmetric and 
nonsupersymmetric one) before vacuum transitions. 

ii) The next main point of the toron approach may be formulated as follows. We 
hope that in the functional integral of the gluodynamics , when the bare charge tends to 
zero and when we are calculating some long range correlation function, only certain field 
configurations ( the toron of all types) are important. In this case the hopeless problem of 
integration over all possible fields is reduced to the problem of summation over classical 
toron configurations. I have no proof that the system of solutions which have been taken 
into account is a complete system. But I would like to stress that a lot of problems ( like 
the 9 dependence , the U(l) problem , the counting of the discrete number of vacuum 
states , the confinement , the nonzero value for the vacuum energy and so on...) can be 
described in a very simple manner from this uniform point of view. 

Both these points are quite nontrivial ones. However , I would like to convince the 
reader in the consistency of these assumptions by considering a more simple models , where 
the answers are well known beforehand. By this reason let me recall in passing that all 
calculations, based on the toron solution [7| demonstrate its very nontrivial role in different 
field theories. Most glaringly these effects appear in supersymmetric variants of a theory. 
In particular, in the supersymmetric CP Ar_1 -theories, the torons (point defects) can ensure 
a nonvanishing value for the < tpip >~ exp(2iirk/N+i9/N) with right ^-dependence. Such 
behavior is in agreement with the value of the Witten index which equals A^jjTJJ and in 
agreement with the large iV-expansion [15]. In analogous way, the chiral condensates can 
be obtained for 4d theories: supersymmetric YM (SYM), supersymmetric QCD (SQCD)( 



see also calculations |16fl,[]17|], based on the standard 't Hooft solution). In these cases a 
lot of various results are known from independent consideration (such as the dependence 
of condensates on parameters m, g; the Konishi anomaly equation and so on...). fl~8| . 



Toron approach is in agreement with these general results. The same approach can be 
used for physically interesting theory of QCD with Nf = N c . In this case an analogous 
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calculation of < ipip > does possible because of cancellation of nonzero modes, like in 
super symmetric theories. For this theory the contribution of the toron configurations to 
the chiral condensate has been calculated and is equal to: < iptp >= — it 2 exp(5/12)2 4 A 3 
[7|. As is well known in any consistent mechanism for chiral breaking a lot of problems, 
such as: the U (l)-problem, the number of discrete vacuum states, the 6*-puzzle, low energy 
theorems and so on, must be solved in an automatic way. We have checked that all these 
properties f?j are consistent with the toron calculation. 

2. The calculation of the topological susceptibility. 

Let me start by giving a few formulae from ref.pl]. The grand partition function is 
given by 



oo \4(kx+k2) fci+fc2 



Z = Y1 71 w; m E II d 4 Xiexp(-e in t.), 



4 2 M 4_1/3 47r 2 

tint. = - o £ <lihln{xi - Xjfqjlj + - lnL 2 {J2 <nhf \ A 4 " 1/3 = c ° exp( — ^tjtt)- 
6 i>j 6 i 9 \ M o) 9 \JV1q) 

where two different kinds of torons classified by the weight Ii of fundamental represen- 
tation of the SU{2) group and qi is the sign of the topological charge. Besides that, in 
formula (|8|) the value g 2 (Mo) is the bare coupling constant and Mo is ultraviolet regular- 
ization, so that eq.(H) depends on the renormalization invariant combination A. As it can 
be seen from (§) the only configurations satisfying the neutrality condition 

E^ = ° (9) 

i 

are essential in thermodynamic limit L — ► oo. In obtaining (|8]) we took into account that 
the classical contribution to Z from k torons is equal to 

4vr 2 

Z ~ exp( —k). (10) 

g 2 

Besides that the factor d A Xi in eq.(§) is due to the 4 translation coordinates accompany 
an each toron f] and combinatorial factor fcilA^! is necessary for avoiding double counting 
for k\ torons and k<i antitorons; lastly, the average overall configurations q, I is an average 
over all isotopical directions and topological charge signs of torons. The constant c in the 
definition of A is the calculable constant 

c = 2Vexp(-a(l)/2) (11) 

where the coefficient a(l) is tabulated in ref. |L9 . 



To compute some vacuum expectation values it is convenient to use the correspondence 
between the grand partition function for the gas (H) and field theory with Sine-Gordon 



5 Let us recall that the one toron has exactly four zero modes in contrary with instanton possessing 
by eight zero modes. 
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interaction , as it was done by Polyakov in ref. [1C] for 3d QED. Let us rewrite @ in the 
form: 

Z e = J D$ex P {- J d 4 xL eff ), □ = d^, (12) 

L eff = 1/2(D0) 2 - A 4 exp(z87r/v / 3/a0 + iO/2) - A 4 exp(-i&ic/V3f a <j> - i9/2). 

In this deriviation it was used the fact that the logarithm function which appears in the 
formula for the interaction (||) is the Green function for the operator □□. After that we 
can use the method [K| to express the generating functional in terms of effective field 
theory (0). 

In this effective field theory the sum over I a runs over the 2 weights of the fundamental 
representation of SU{2) group. Note, that the first interaction term is related to torons 
and the second one to antitorons. Besides that, since we wish to discuss the 9 dependence 
, we also include a term proportional to the topological charge density -^^G^G^ to the 
starting lagrangian and corresponding track from this to the effective lagrangian(|T2"D. 

The most important result from ref. || is the nontrivial dependence on 9 of the topo- 
logical density and susceptibility.These quantaties are relevant for the solution of the U{1) 
problem : 

(^G^G~) = iF{9) = i2K A sin(0/2), -tt < 6 < n. (13) 



*327T 



J d 4 xexp{ikx)(^-^G f , u G f , u (x), ^-^G^G^O))^ ~ ~ ~ cos(^). (14) 

As was discussed in ref.|l2|, the reason to have the nontrivial 6- dependence (|13|'0) as 
well as the strong quasiparticle interaction ~ ln(xj — Xj) 2 (§) in YM theory, is the presence 
of the nontrivial algebraic structure ~ q^iq^j in the expression for e int (H). Just this 
fact was crucial in the analysis of the 2 dimensional CP^ -1 model also [Q.Such structure 
for ei n t is in the striking contrast with 2+1 Polyakov's model [f[0|| , where the interaction 
energy proportional to the topological (magnetic) charges of quasiparticles only: 

tint ~ | r- (15) 

Just this difference leads to the existence of the nontrivial 9 dependence in 4d YM theory 
(in spite of the fact of the strong quasiparticle interaction ~ ln(xj — Xj) 2 ) in contrary with 
Polyakov's model where physics does not depend on 9. 

We proceed now to the direct calculation of the topological susceptibility ([14]). 

From the eq.(^) it is easy to reproduce the formula ( |i~3"D for the vev of the topological 
density. To this aim we differentiate the Zg with respect to 9 and put <fi = in the 
vacuum state. To compute a correlation functions depending on x 2 we have to introduce 
2 different fields Xai & — 1> 2 in place of the 4>: 

Xa = ~^ I » ( t > - ( 16 ) 
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These Xa fields should be considered as independent in the course of calculation , because 
the constraint 

J2Xa = 0(mod2n) (17) 

should be set to zero only at the end of calculation. The reason for this is as follows. 
The constraint (|T?f ) is appeared as a manifestation of the neutrality condition (0)in the 



infinite volume limit of the functional description (|12|) for our ensemble (|8|). We should 
be very careful in taking such limit in the course of calculation of a different topological 
characteristics. It can be argued that this limit should be set at the last stage of calculation 
only ( otherwise, we would get some senseless results). Thus , the constraints ( |17"D as well 
as (||) will be fixed at the end of calculation. 

With the above consideration taken into account , the topological susceptibility can 
be expressed in terms of the new variable Xa'- 

f 1 ~ 1 

Xt{k) = d 4 xexp(ikx)(—^G^ u G^ u (x),——rG tll ,G llu (0))= (18) 



-A./^«p(*x)<-nxT»),-nX I (0) + ^(x,,-nx,(0)> 

Thus, the problem of calculation of the topological susceptibility reduces to the problem 
of calculation of the following correlation function: 

(M+( Xl ),M a (x 2 )), M a (x)=exp(t Xa ). (19) 

where the operator M is so called disorder operator. 

The importance of disorder variables in gauge theories has been emphasized by 't 



Hooft p0 |, who has argued that rather than instantons it is the field configurations with 
nontrivial Z^ topological charge that should be considered responsible for the long range 
confinement of quarks. Analogous disorder variables have been used in different fields 
of physics [0J.In the context under consideration this variable was introduced in 2+1 
Polyakov's model in ref . [ 22| and in two dimensional CP N ~ l model in ref.H. 



A set of operators is now defined as follows. M is an operator that acts on original 
fields by gauge transforming them by U x °(x); this gauge transformation is singular at 
Xo and has the property that as x encircles xq,U does not return to its original value ( as 
it happens in the instanton case), but acquires a phase (N = 2 for SU{2) group): 

U X0 ((f) = 2tt) = exp(-i2n/N)U X0 ((j) = 0), U = exp(2iT0) (20) 

where <fi is an angle variable in the x — y plane and the point xq lies at the same plane. 
From its definition it must be clear that M(x) absorbs one half topological unit, so we 
say that M(x) is the annihilation operator for one point toron at x with weight I and 
M + (x) is the creation operator for one toron.Q It should be clear , that U depends on all 

6 Let me note that the small size instanton (p — > 0) at x$ is also can be created by the singular gauge 
transormation U(x). However, in the instanton case U does return to its original value ( in contrast with 
(p0|)) as x encircles xq in the xy plane. 
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x M variables , so that M is the annihilation operator for the point defect. However, at the 
x — y plane, U depends only on angle variable <fi. The singularity of = iU + d fl U must 
be smeared over an infinitesimal region around x as it was done for the separate toron 
solution^]. We will not consider the regularization problem in this paper. 

Now we would like to demonstrate that the definition of M as an operator of large 
gauge transformation in terms of original fields exactly leads to the expression (|l9|)in 
terms of the effective field theory ([12]). To this aim, let us consider the e, n £ in the formula 
(|) after the action of the gauge transformation (|20|) at point x . Because this gauge 
transformation creates an additional toron at point xq with isospin Iq in the system of 
the other torons placed at X{ with isospins l{ we will obtain an additional contribution to 
the €i n t- Namely, after action of the operator M we have an additional interaction term 
between created toron Iq and torons J, from the system 

Ae int ~ fofi M^o - Xi) 2 . (21) 

i 

It is easy to understand that this interaction after simply repeating the deriviation of 
eq.(|12D, reduces to the following expression in the effective field theory: 

(M(x )) = J D$exp{- J d 4 xL e// .) exp(z87r/V3/^(^ ) +W/2), □ = 0^,(22) 

L eff = 1/2(D0) 2 - ^ A exp(i8n/VsfJ+i9/2) -Y / ^ 4 exp(-t8n/V3fJ-i9/2). 
f-+± f 

Thus, the operator M under consideration in the effective theory looks like that 

M a (x) = exp(z87r/V340(i) + iO/2) = exp(i Xa {x) + id/2), (23) 



where we rewrite eq. ( P3|) in terms of the x fields (0) . 

I would like to emphasize , that the expression for the operator of the large gauge 
transformation in terms of the effective variables (|23| ) has the same form like in another 
field theories, 2+1 dimensional QED |22| and 2 dimensional CP N ~ l models H. As it was 
argued in the last reference, it is not an accidental coincidence, but this fact has a very 
general origin. 

Now we are in position to calculate the correlation function 

d 4 xexp(ik(x 1 -x 2 ))(M(x 1 ),M + (x 2 )) (24) 



As it was explained in refs.[ID], |22| in the analogous calculations in 3 dimensional QED 
we can not neglect terms related to insertions of M(x). Rather, this functional integral 
can be estimated by the steepest descent method. Thus, if we assume that this functional 
integral is dominated by the classical field, \ satisfies 

Att Atx 

OO x + 2A 4 (- ? =) 2 sin X = ~ xi) ~ - x 2 )). (25) 



7 Usually we have in mind the special kind of regularization which preserves the self-duality equation, 

Q 2 

so that the classical action for this configuration equals 
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The physics suggests linearization is legitimate so we can substitute x instead of sin(x) in 
the equation (|25| ) and find Xci{ x ) by means of Fourier transformation. However, for our 
purpose an explicit form of Xci(%) is not needed , because we are interesting in the Fourier 
transformed form of the correlator ([24]). Simple calculation gives the following expression 
for this correlation function at k — ► 0: 

<fxe W {ik{x 1 - x 2 ))(M(x 1 ), M + {x 2 )) ~ (A: 4 -^ + A 4 ) -1 . (26) 

Now several comments are in order. First, in comparision with the well known expres- 



sion for topological susceptibility at k = (II ) we were able to find the k dependence at 



small k. It turns out that the k dependence comes in this formula through k 4 . Such be- 
havior, from the one hand, is the direct consequence of the strong logarithmic interaction 
of the pseudoparticles (||) and from the other hand it implies that 

x ' t (k = 0) = A J d*xeMikx)(^G^G^(x), -L^G^G ^(0)) k=0 - 0. (27) 

If such behavior were in QCD theory , this would mean the zero magnitude for the singlet 
axial form factor, Ga{0) [l3j m agreement with EMC experiment. However, this is not 



the case, because we have dealt with the pure YM theory without quarks and one could 
expect that the dynamical light quarks will change this situation. But we may expect that 
such k A behavior can not disappears without leaving a trace in the full theory. However , 
we can only speculate on this topic now. 

As a second remark, we note that the link between the 9 dependence in the the- 
ory, discrete number of vacuum states there and confinement phenomenon mentioned in 
Introduction can be understood by the following way. 

From the one hand, we have nonzero value for the vev of the disorder operator M (|23|) 

(M) ~ exp(2i7rk/N + i6/N), = 0, 1, 2..., N — 1. (28) 

In obtaining eq.(^) we took into account that in different vacua the fields Xa take N 
different vacuum values Xa = 2nk/N , as was discussed in ref.||0. On the classical 
level we have Z N degeneracy , but on the quantum level, taking into account the toron 
vacuum transitions, these vacua have different energy as was the case in the standard 
\9 >= yjexp(in6 l )|ri > vacuum consideration, but now there is only a finite number of 
vacua.The operator M plays the role analogous to the chiral condensate in QCD. Different 
worlds discussed above are labeled by the phase of this disorder parameter M. 

From the other hand, one could expect that < M >^ just corresponds to confine- 
ment phase. I have no rigorous proof for this, however , an analogous considerations in 
2 dimensional CP^N — 1) model M and in 3 dimensional Polyakov's model [22|, where 



vev of the disorder operator have been calculated ,verify this conjecture. Moreover, the 



8 It is easy to understand that this is the direct consequence of the periodicity of the effective lagrangian 
(|l2l). Although each of the separate solutions has a 9 period of 27riV,the overall minimum has a period 2n 
because of jumping from one value to another at 9 = ir. 
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explicit calculation of the vev of the Wilson line in 4 dimensional gluodynamics (see the 
next section), also demonstrates the nontrivial relation between nonzero value for < M > 
and area law for Wilson loop < W >. 

As a last remark, I would like to note, that the correlator (Jl8|), describing a gluonium 
states, can be expressed in terms of the operator M. The corresponding correlation 
function in the x space 

(Af(xi), M + (x 2 )) (xi . X2) ^ QO - (M) 2 ~ exp -\ Xl - x 2 \ (29) 

shows an exponential fall off at large distances. Thus, one could expect that the oper- 
ator M can be identified with the creation operator of the lowest gluonium state. It is 
interesting to note that a solitonic operator in the different field theories as well as the 
vertex operator in the conformal field theories have the same exponential form. Thus, 
the creation operator M which should be highly nonlocal and nontrivial in terms of the 
original fields (A^-gluons) has a very simple form in the terms of the auxiliary (dual) 
variables Xa- 

3. Wilson loop operator 

In order to show directly ( without references to the analogy with 2 and 3 dimensional 
systems) that there are electric strings in the ensemble of our pseudoparticles @ let us 
calculate the vev of the Wilson loop. As was explained above, the torons at large distances 
look like singular pure gauge field with definite isotopical direction and so, the A® field 
is abelian at large distances (in a more detail see ||). Thus, the standard quasiclassical 
approximation, when we substitute for the corresponding classical solution, leads ( 
after simply repeating the deriviation of (0)) to the following expression for < W > at 
# = 



(W) = {TrexpifiqApdxJ) = (30) 
J D<pexp(- J d 4 xL eff .), L eff = 1/2(D0) 2 - £ 2A 4 cos(8tt/V3X<? + f a $). 

where term proportional to $ is related to Wilson loop insertion and has the following 
property : $(x) is equal to the external charge 27rg if x e S, Wilson plane, and $ = 
otherwise. In this deriviation we took into account that if the toron is in the S plane, 
then the integral over G^da^ is non-zero, and it is equal to zero otherwise (see for a 
more detail about toron properties the ref .0 ). 

Since the field $(x) is strong enough , we can not neglect non-linearities in ( |3~0D just like 
in the calculation of the topological susceptibility (see previous section). The estimation 
of this functional integral can be done as before by the steepest descent method with 
respect to <fi. The corresponding field <p c i is determined from the classical field equation, 

4-7T 

nn x ' + 2A 4 (^) 2 sin(x') = 2k0s(z, t)S'(x)5\y), (31) 

where \' = X + ^I-> 5'(x) = ,9s(z,t) — 1 if z, t £ S and 9g(z,t) = otherwise. The 
right- hand side of this equation is related to the Wilson loop insertion, i.e. with function 
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$(x). Indeed, for the Wilson loop placed in the t — z plane we have $ ~ 5 Xi o^y,o an d thus 
, acting by the operator □□ on $ we exactly reproduce the right hand side of the eq.fl3T|) 
( without loss of generality we consider x > 0,y > 0): 

□□4,0^,0 ~ □□(! - e(a:))(l - e(y)) ~ n[5'(x)(l - e(y)) + 5'{y){l - e(x))}, (32) 

<9 2 d 2 

where we have substituted the expression 1 — e(x) in place of S X) o and used the following 
feature of the e(x) function :Oe(x) ~ S'(x). The next step of calculation of the r.h.s. (|32"D 
is the using of the Fourier representation for the functions e(x) ~ J °° sm ^ x dk and S(x) ~ 
J^° oa exp(ikx)dk which leads to the desired result: 

r.h.s.(32) ~ □ f dk x k x [l - - f dk y —^^- \ exp(ik x x) + (x -> y) (33) 



dk x k x exp(ik x x) / dk y k y sm(k y y) + (x — > y) ~ 8'{x)8'{y). 

oo JO 

Now we are ready to estimate the functional integral (|3CiD by the steepest descent 



method. For this purpose , we substitute x instead of sinx in the eq.(|3TD and find Xd 
by means of Fourier transformation just as it was done before in the previous section 
for the calculation of the topological susceptibility We expect that the accuracy for such 
procedure is not very high and the numerical coefficient given bellow should be considered 
as an estimation of the order of value. With these remarks in mind we obtain 



(W) ~ exp(-^A 2 S) (34) 
bArc 



The result Q53| ) implies that between two fixed charges there exists an electric string with 
energy density ~ A 2 . 

Now several comments are in order. First of all, I have to stress that the unusual kinetic 
term (Ox) 2 in the effective action (|TJ)plays a crucial role in this deriviation. From the one 
hand ,such kinetic term is the direct consequence of the strong pseudoparticle interaction 
(§) ,and from the other hand just this term gives a correct expression ~ 5'(x)5'(y) for the 
right hand side of eq. (|3l"l) , leading to the existence of the solitonic shape solution in the x, y 
directions. ^ This solitonic shape guarantee the convergence of the corresponding integral 
/ dxdy over x, y directions, perpendicular to Wilson plane z, t and leads to confinement 
formula (pip. The thickness of the sheet with Wilson loop as boundary is of order A -1 , 
plasma screening length, and not A — > 0, toron size. This is just the reason to get the 
finite answer for the string tension ( p4|) in comparision with our very rough previous 

9 Let us note that the analogous problem in 3 dimensional QED looks as follow . If Wilson loop is placed 
in the xy plane, the equation , analogous to ( |3T| ) is essentially one dimensional ]To| and (xc/)depends only 
on z . The solitonic shape for the corresponding solution is due to the non-linearity (~ $'{z)) related to 
the Wilson loop insertion. This solitonic shape guarantee the convergence of the corresponding integral 
J dz, over perpendicular to Wilson loop direction and leads to the confinement formula. 
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estimation ||, where back reaction of the Wilson loop insertion has not been taken into 
account. 

As a second remark ,we note that our explicit calculation of the < W >~ exp(— S) 
is in the consistency with our conjecture that the confinement phase just corresponds to 
the nonzero value of disorder operator < M > fl28|) . As it was mentioned above, it is 
also correct in the 2 dimensional CP 1 ^^ 1 model || and 3 dimensional QED [22| where an 
explicit calculations do possible. 

4. Final remarks. 

The main point of this Letter is that the dynamical solution of the one from the 
following problems : U(l) problem, confinement problem, multiplicity of the vacuum states 
N, 6/N dependence, should be necessarily accompanied by the resolution of the rest 
problems within the same approach. We have checked this conjecture in the framework 
of the dynamical toron calculation. I expect that two main assumptions of this approach 
I am using 

i) the multivalued functions are admissible in the definition of the functional integral; 

ii) the only certain field configurations ( torons) are important and the problem of 
integration over all possible fields is reduced to the problem of summation over classical 
toron configurations, are consistent assumptions. I have no proof for this, but I would like 
to stress that all problems mentioned above can be considered from this uniform point of 
view. 

The main quantitative results, which have been obtained can be formulated as follow: 

^-OW^+AV (35) 

< M >~ exp(z6>/2 + iirk),k = 0, 1. (36) 
< W >~ exp(-aS) (37) 



The first eq.(|35|) shows that the k dependence of the topological susceptibility comes 
through k 4 . Such behavior might be important for explanation of EMC experiment. 

The second eq.(|36|) explicitly demonstrates that vacua in YM theory are classified by 
the phase of the disorder operator M. Apart from standard parameter 6, < 8 < 2ir, 
there is a new label k = 0, 1. Although each of the separate solutions has a 6 period 
47r, the overall solution has a 9 period of 2ir because of the jumping from the one value 
(k = 0) to another (k = 1) at 9 = n. Such behavior is in agreement with large N results, 
where the number of vacuum states is equal N and 9 dependence appears in the form 
6/N. Besides that,< M >^ is in consistency with the confinement property of the 
theory. 

The explicit calculation of the vev of the Wilson line ( |3~7|) also confirms this conjecture. 
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